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NAVIGATION DESIGN AND ANALYSIS FOR THE ORION 
EARTH-MOON MISSION 

Christopher D’Souza* and Renato Zanetti f 


This paper details the design of the cislunar optical navigation system being pro- 
posed for the Orion Earth-Moon (EM) missions. In particular, it presents the 
mathematics of the navigation filter. The unmodeled accelerations and their char- 
acterization are detailed. It also presents the analysis that has been performed to 
understand the performance of the proposed system, with particular attention paid 
to entry flight path angle constraints and the AV performance. 

INTRODUCTION 

Vehicles navigating to or from the Moon usually rely on ground tracking and ground updates to 
perform the insertion and correction maneuvers. A natural advancement in technology is autonomy. 
The Orion vehicle, designed to explore space beyond LEO, is required to return the crew safely 
in the case of loss of communication with the ground. As such, it needs to be able to navigate 
autonomously, independent of ground-based measurements, utilizing on-board sensors subject to 
stringent mass/power/volume constraints. Since the vehicle will be carrying optical cameras, the 
cislunar navigation system is designed to use images obtained from these cameras, in particular 
star/planetary limb and planetary disk measurements. Whereas the navigation system of Orion in 
and below LEO is well understood, the design of the cislunar navigation system unique presents 
challenges. 

Whereas the Orion sensor complement includes two star trackers, the star trackers being con- 
sidered have a very limited field-of-view. As such, they don’t lend themselves to cislunar optical 
navigation, which needs fields-of-view in excess of 20 degrees. Thus, optical cameras, which arc 
already planned for situational awareness, arc harnessed into a cislunar navigation role. In this pa- 
per the design of the cislunar optical navigation system being proposed for the Orion Earth-Moon 
(EM) missions is presented. In particular, it will present the mathematics of the navigation filter and 
the analysis that has been performed to understand the performance of the proposed system, with 
particular attention paid to entry flight path angle constraints and the DV performance. 

Previus studies focused on the lunar orbit determination problem [1,2]. Tuckness and Young 
consider autonomous navigation for lunar transfers [3]. Their analysis focuses on azimuth and ele- 
vation measurements of the Earth. Moon, and Sun. Two star-elevation measurements relative to the 
planet’s limb provide the same kind of information as azimuth and elevation of the apparent center 
of the planet. However, star-elevation measurements arc the preferred approach for two reasons. 
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First, multiple stars can be processed simultaneously, and the redundant information effectively fil- 
ters out some noise. Second, the method used here does not depend on the attitude of the spacecraft, 
or on the misalignments of the sensors, reducing the possible error sources. 

This work builds on previous results presented in [4] but focuses on new analysis done to char- 
acterize the process noise, i.e. unmodeled accelerations, and the preliminary results of the EMI 
navigation system design. Christian and Lightsey ADD CITATION AND DESCRIPTION 

This investigation specifically addresses the transfer from the Moon to the Earth. In an emergency 
situation, during a loss of communication scenario, the primary objective is the safety of the crew. 
This subsequently translates into a flight-path angle requirement at entry interface (El) for a direct 
entry. A direct entry, as opposed to a skip entry, reduces the risk of the capsule bouncing back into 
space, and allows for a greater margin on the flight-path angle at El. 

The accuracy of the flight-path angle at El is driven by several factors including the navigation, 
targeting, and burn execution errors at the time of the last mid-course maneuver, and unaccounted 
trajectory perturbations between the last mid-course maneuver and El. Apollo missions tolerated a 
maximum flight path angle error at El of ±1 degree, with half of this error allocated to navigation. 
A similar criterion is employed in this study. 

Perturbations are a major source of errors in the cislunar navigation performance of Orion. In a 
perfect world all the sources of perturbations would be modeled in the filter dynamics. However, 
computational limitations preclude such extensive modeling. Therefore, the primary sources of 
perturbations are characterized. In particular there are three categories of unmodeled acceleration: 
propulsive sources, gravitational perturbations, and solar radiation pressure. Only propulsive errors 
are included in this analysis; the gravitational and solar radiation pressure are not included - they 
will be included in a future study. For EMI, the gravitational and solar radiation pressure errors are 
several orders of magnitude below the thrusting sources. The propulsive sources considered are: 
attitude deadbands, attitude slews, C02 venting, and sublimator venting. 

Lineal - covariance techniques are used to perform the analysis for the Orion Cislunar missions. 
This comports well for the navigation system design since the cislunar navigation system on Orion 
will be an Extended Kalman Filter. Many of the same states and dynamics used in the linear co- 
variance analysis will be used in the on-board cislunar navigation system. A preliminary design of 
the cislunar navigation system is presented. This is supported by linear covariance analyses which 
provides navigation performance, trajectory dispersion performance and AV usage. 

The paper is organized as follows: Section 2 will contain a brief description of linear covariance 
analysis. In Section 3, the navigation system will be described. Section 4 will contain a description 
of the perturbations used in this analysis. Section 5 will contain results of this analysis. Finally, 
some concluding comments will be made in Section 6. 


LINEAR COVARIANCE ANALYSIS 


This investigation is performed using linear covariance (LinCov) analysis techniques [5,6]. The 
state vector is given by [4] 


x 


= I 1 
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The nominal trajectory is obtained by integrating the nominal dynamics model with an Encke- 
Nystrom method [7]. Neither the rotation vector 6 nor its uncertainty are integrated in this analysis. 
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The nominal attitude is known at any time and it does not need to be calculated. The attitude 
estimation error covariance is constant and is driven by the star tracker accuracy. The attitude 
navigation dispersion covariance is constant and is given by the attitude control dead-hand. The 
attitude environment dispersion covariance is constant and obtained from the above two quantities 
assumed uncorrelated. Before the star elevation is determined, the vehicle slews in preparation 
for measurement acquisition. This attitude maneuver is performed by the onboard thrusters and is 
assumed to be instantaneous. Due to thruster misalignment, this maneuver adds uncertainty to the 
translational states. After the batch of measurements is available, the vehicle returns to its nominal 
attitude. In linear covariance analysis, the difference between the true state and the nominal state is 
defined as the environment dispersion 

fix = x — x. (2) 

The difference between the estimated state and the nominal state is defined as the navigation disper- 
sion 

fix = x — x. (3) 

Finally, the difference between the true state and the estimated state, is defined as the estimation 
error, sometimes referred to as the onboard navigation error 

e = x — x. (4) 

Following the standard Kalman filter assumptions, the difference between the nominal and estimated 
models is represented with zero-mean, white noise. The estimated state evolves as 


x = f(x), (5) 

where f is a nonlinear function representing the system dynamics as modeled by the filter. The 
evolution of the nominal state is modeled as 

X = f(x) = f(x) + v , (6) 

where f is a nonlinear function representing the state dynamics as modeled in designing the nominal 
trajectory. The nominal dynamics f may be higher fidelity than the filter’s dynamics f. The vector 
v represents the dynamics modeled in the nominal trajectory but neglected in the filter models. 
In Kalman filtering, the difference between the true dynamics and the filter’s dynamics is called 
process noise. While these unmodeled dynamics are not actually white noise, they are modeled as 
such. The power spectral density of process noise is then tuned to achieve good performance. The 
same procedure is used here. In order to capture the difference between the two dynamical models, 
v is modeled as a zero-mean white process with spectral density Q. The goal is to represent the 
increased value of the navigation dispersion during propagation due to the difference between the 
nominal and filter’s dynamical models. 

The evolution of the navigation dispersion can be approximated to first-order as 

fix = x — x = f(x + fix) — f(x) — v ~ F(x)5x — v. (7) 

The evolution of the navigation dispersion covariance is governed by 

P = F(x)P + PF(x) t + Q. (8) 
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Similarly, the true state is modeled to evolve as 


X = f (x) + V. (9) 

The evolution of the estimation error is given by 

e = x — x ~ f(x) + F(x)(x — x) + v — f(x) — F(x)(x — x) = F(x)e + v. (10) 

Vector r/ is modeled as zero mean white noise with spectral density Q. The onboard covariance P 
evolves as 

P = F(x)P + PF(x) + Q. (1 1) 

Notice that the Jacobian F could be evaluated at the estimated state x instead of the nominal state 
x, as in the extended Kalman filter. 

Finally 

5x = x — x ~ F(x)<5x + v — v (12) 

and P evolves as 

P = F(x)P + PF(x) + Q. (13) 

Notice that Q = Q + Q if u and v are assumed to be uncorrelated. 

Since the environment and navigation dispersions arc naturally correlated, it is intuitive to create 
an augmented dispersion state, whose covariance is defined as II 
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The evolution of the augmented covariance is given by 
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where it is assumed that u and v are uncorrelated. All error states are modeled as first order Markov 
processes and arc assumed to be uncorrelated to each other. 

THE ORION CISLUNAR NAVIGATION FILTER 
The Filter Dynamics 

Since this filter operates once Orion is away from Earth (and outside of GPS range), the primary 
forces governing the motion of the vehicle arc the gravitational forces of the Earth, the Moon, 
and the Sun. The trajectory is designed taking into account all three of these bodies. Whereas 
the equations of motion arc formulated with respect to a central body, this (central body) changes 
depending on which sphere of influence the vehicle is subject to. 

The equations of motion for the Earth-Sun-Moon system are 
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where and r py is the position of the vehicle ( V) with respect to the primary body (P), r qv is the 
position of the vehicle with respect to the secondary body ( Q ), r pq is the position of the secondary 
body with respect to the primary body, rgy is the position of the vehicle with respect to the Sun (S), 
and rps is the position of the Sun with respect to the primary body. The geometry is shown in Figure 
2. In many applications, these equations are integrated by a Runge-Kutta or Runge-Kutta-Fehlberg 
fixed-step or variable-step algorithms. 

However, blindly applying a standard fourth-order Runge-Kutta method can lead to numerical 
errors (if large step sizes are taken) not to mention inefficiencies (if small step sizes are taken). The 
Encke-Nystrom method has been shown to have none of these deficiencies - large step sizes can be 
taken and numerical precision maintained at the same time. 



Figure 1. The Earth-Sun-Moon- Vehicle Geometry 


The Encke-Nystrom Method The Encke-Nystrom method has an illustrious legacy. Unlike the 
Cowell’s method, only the perturbations away from two-body motion are integrated. As such, the 
perturbations being integrated arc small and excellent numerical precision is retained. In addition, 
since the force-field under which the equations of motion being integrated is conservative, and 
since the velocity derivatives are expressed in terms of the position only, a Nystrom formulation of 
integration is used. The disadvantage of this method is that there is a bit of mathematical set-up 
required. However, if one is willing to pay this price, the dividends arc enormous - both in terms of 
precision and computation time. The equations which are integrated are as follows: 

tpv = -"if-rpv ~ [ /(qq)^pq + rpy] - [f(qs)*ps + r PV ] (18) 

r PV r QV r sv 
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This formulation was used in the Apollo program [7], It should be noted that Eqs.(17) and (18) arc 
mathematically equivalent. 

The Trajectory Partials In order to compute the state transition matrix from an epoch to an- 
other epoch in the future (or the past), the partial of the future (or past) state with respect to the 
initial epoch must be computed. The reason the Encke-Nystrom algorithm is introduced into this 
discussion is that the Encke-Nystrom method is predicated on “small” deviations from a reference 
(osculating) orbit. This is precisely the foundation upon which the state transition matrix for the 
cislunar trajectory is built and it is also the foundation upon which linear covariance analysis is 
built. Whereas the state transition matrix is defined as 
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With the dynamics defined in Eq. (40), the partials of the dynamics are found to be 
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where G(t) is defined as 
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The Measurements 


Optical measurements were processed every 60 seconds. The optical camera had a 18 degree 
Field of View. The optical measurements consist of star-horizon measurements and apparent angular 
radius measurements. These measurements take into account when the planet is both larger than and 
smaller than the field of view. 

UNMODELED ACCELERATION CHARACTERIZATIONS 

There arc three categories of unmodeled acceleration with which a vehicle in cislunar space 
has to contend: propulsive sources, gravitational perturbations, and solar radiation pressure. Only 
propulsive errors arc included in this analysis; the gravitational and solar radiation pressure arc 
not included - they are several orders of magnitude below the propulsive sources. The propulsive 
sources, therefore, can be further separated into thruster errors and ECLSS (Environmental Con- 
trol and Life Support System). Thruster errors include attitude deadbanding and attitude slewing 
maneuvers. ECLSS sources include Pressure Swing Adsorption (PSA) ( CO 2 venting), ammonia 
sublimator venting, and waste water venting events. EMI will have a ‘metabolic simulator’ which 
will exercise the entire ECLSS system except for waste water vents; hence waste water venting 
perturbations arc not included in this analysis. 

Attitude Dead-banding Maneuver Errors 

Given a 24-jet ESA-SM config, with attitude dead-hands with 0.028 second on-time, the expected 
ALis 

' —8.2E — 4 ' 

AV = 8.2E — 4 ft/sec (20) 

1.0E — 20 _ 

and assume that there is one firing every 30 minutes*. This means that there is a la AV (assuming 
no inertial direction preference) of 

a a v,rss = 1-160 x 10~ 3 ft/sec a av, p er axis = 6.695 x 10~ 4 ft/sec (21) 

This gives a process noise of 

QAtt Deadband = = 2 ‘ 490 X 10 ~ 10 ft VseC 3 (22) 

1800 sec 

Attitude Slew Maneuver Errors 

For attitude maneuvers, assume 25 attitude events, resulting in 50 slews over the mission, given 
a mission time of 160 hours. This can be thought of (and spread out) as one attitude slew every 
3.2 hours. These events are as follows: 7 maneuvers (14 slews), 7 optical passes (14 slews), 10 
communication passes (20 slews), 1 other attitude event (2 slews). 

Given a slew maneuver with a 0.7 second on-time with a 24 jet configuration, the residual error 
in body coordinates is 

' —2AE - 2 ' 

AV = 2AE-2 ft/sec (23) 

L0 E-A 

"Rodolfo Gonzales, 4/12/2013 
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The 1 sigma error is 


&AV,RSS, Slews = 3.394 X 10 2 ft/sec 4=>- CF AV, Per axis, Slews = 1-960 X 10 2 ft/sec (24) 

So, the process noise for attitude events with a 24-jet configuration is 

2 

-V \ P er axis , Slew s n oor- , , a r \ — 8 r. — / 3 

Q AttSlews = 11cnn = 3.335 x 10 ° ft /sec (25) 

11520 sec 

For the case of continuous tracking, with a pan/tilt camera, there will be less need for attitude 
slews to perform celestial navigation. Hence it is assumed that there will be 7 fewer attitude events, 
for a total of 18 attitude events, resulting in 36 slews over the mission, given the same mission time 
as before, resulting in one slew every 4.44 hours. The process noise is 

Qausi = a ' iv : P ;Z‘‘' S,e '‘‘ = 2 - 401 X lcr8 ft2 /sec 3 (26) 

16000 sec 

For the case of 12 attitude events (24 slews), this will result in one attitude slew ever 6.66 hours. 
So, the process noise for attitude events with a 24-jet configuration is 

Qaus ira , = = 1.601 X 10- 8 fV/sec 3 (27) 

24000 sec 

Finally for the limiting case of only 6 attitude events (12 slews), we will have one attitude slew 
every 13.33 hours, resulting in a process noise of 

Qausi™, = = 8-003 X lcr* ftVsec 3 (28) 

48000 sec 


PSA Puffs 

The Pressure Swing Adsorption (PSA) cycles automatically every 6-10 minutes, less during crew 
sleep periods, more when the crew is awake. The PSA system can be commanded to force a com- 
plete desaturation (all vents open into space) and then they would likely not open automatically for 
40-60 minutest This would be needed during ground tracking passes and optical navigation passes. 

The AV resulting from the PSA puffs can be calculated as follows: the force during a single vent 
is 36 lb-force and it lasts for 0.1 sec with a spacecraft mass of 1650 slugs. Therefore, the AV is 

FAt , , 

AVpsa = = 2.182 x 10 3 ft/sec AVpsAneraxis = 1-256 x 10 3 ft/sec (29) 

m 

Given that this will run every 6 minutes, 

2 

x, ® AV, PS A, Per Axis . Ar , a , n _9 r , 9 / 3 /on\ 

Qpsa = VpF = 4.408 x 10 fG/sec 3 * (30) 

360 sec 


^Email correspondence from John Lewis and Chris Delnero 



The Quiescent Process Noise 


The combination of these three events is 

Q^Quiescent = QpSA + Q%fsiews + Q Att Deadband = 3.801 X 10 8 ft 2 /sec 3 (31) 

This results in 

QQutTcent = 3-531 x 10~ 9 m 2 /sec 3 (32) 

For the case of a pan/tilt optical camera and 36 slews, the quiescent process noise is 

Q Quiescent = QPSA + QAttSUws + Q AttDeadband = 2-8667 X 10~ 8 ft 2 /sec 3 (33) 

This results in 

Q 3 Qu£7cent = 2-663 x 10~ 9 m 2 /sec 3 (34) 

For the case of a pan/tilt optical camera and 24 slews, the quiescent process noise is 

QQutTcent = QpSA + (& + Q AttDeadband = 2.067 X 10" 8 ft 2 /sec 3 (35) 

This results in 

QQutTcent = 1-920 x 10~ 9 m 2 /sec 3 (36) 

For the case of a pan/tilt optical camera and 12 slews, the quiescent process noise is 

QQutZnt = QPSA + QAUsTels + Q AttDeadband = 1-266 X 10“ 8 ft 2 /sec 3 (37) 

This results in 

QQutTcent = 1-176 x 10~ 9 m 2 /sec 3 (38) 


The Active Process Noise 


Finally, there is the case of the ammonia sublimator, which runs occasionally. The sublimator 
vents ammonia ( NH 3 ) when the vehicle is close to the moon (due to the increased thermal environ- 
ment) and after CM/SM separation. When it operates, it imparts a AV on the vehicle. It is expected 
that the sublimators will run for 1/2 hour behind the moon and 1/2 hour before El. For the CM/SM 
configuration. 


AV = 


F At 0.81bf x 1800sec 


m 


1650slugs 

For the CM only (entry) configuration, 

. FAt 0.81bf x 1800sec 
AC = = 


m 


673slugs 


0.873 ft/s = 0.266 m/s 
2.140 ft/s = 0.652 m/s 


This means that 


Q sublime, 


CM/SM 


Q sublimec M 


only 


{ 0 266 \ ^ 

/1800sec = 1.310 x 10 -5 m 2 /sec 3 

( 0.652\ 2 

— /1800sec = 7.877 x 10 ~ 5 m 2 /sec 3 


(39) 


(40) 


(41) 

(42) 


It may well be that we will be able to sense the sublimator vents; if this is the case, then the process 
noise would be reduced, though this would be the case only for the on-board navigation state; the 
ground-based navigation solution would not be affected. 
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Summary of Unmodeled Accelerations 

Table 1 contains a summary of the previously described unmodeled acceleration as well as their 
relative strengths. 


Table 1. Type and Strength of Unmodeled Acceleration 


Type of Noise 

Assumptions 

Strength ( m 2 / s 3 ) 

Attitude Deadbanding 
Attitude Slewing 
Attitude Slewing 
PSA Vents 

Ammonia Sublimator 
Ammonia Sublimator 

Jet firing every 30 minutes 
50 attitude events 
25 attitude events 
Every 6-10 minutes 
In Lunar Vicinity (0.5 hour) 
0.5 hour prior to El 

2.313 xlO -11 
3.098 xlO -9 
1.601 xlO -9 
4.095 xlO" 10 
CM/SM: 1.310xl0~ 5 
CM Only: 7.877 xlO" 5 


RESULTS 

The trajectory under consideration a 8 day trajectory departing on December 15, 2017. It is shown 
in Figure 2. 

For this analysis, we included three sets of navigation concepts of operations: 1 hour optical 
passes, continuous optical passes, and ground tracking only passes. The second case (continuous 
optical navigation) serves as a bounding case to demonstrate the best performance of the navigation 
system were we able to get continuous measurements. Currently, it is not clear where and how the 
optical camera will be mounted. If it is not mounted on a pan/tilt platform, there will need to be 
attitude maneuvers to obtain optical imagery for navigation. For the first case, we assumed that 
there would be one hour of optical passes to the closest body; due to thermal and other operational 
constraints, we will not be able to point the cameras at the nearest body for more than an hour. This 
is expected to be the most reasonable optical navigation concept of operations for this scenario. In 
this case, we would commence tracking two hours prior to each maneuver and terminate the pass 
one hour prior to each maneuver. The concept of operations for ground tracking is as follows: after 
a period of tracking by ground stations, a navigation state (with associated covariance matrix) will 
be uplinked to the spacecraft one hour prior to the maneuver. This state would be used to compete 
and perform the maneuver. This computation of the maneuver could either be done on-board or on 
the ground. 

The optical camera was assumed to have a noise of 5 arc-seconds (ler) and a bias of 3.33 arc- 
seconds. In addition the stellar sub-point had a noise of 10 arc-seconds (ler) and a bias of 5 arc- 
seconds. Finally the horizon had a noise of 10 km (ler) and a bias of 3 km. 

We assumed that there were 7 Trajectory Correction Maneuvers performed. They were chosen 
reflective of what was done during Apollo and will likely be adjusted, but they serve as reasonable 
place-holders, taking into account factors such as crew sleep and time since the previous maneuver. 
In order to ensure accurate delivery a maneuver was placed 6 hours prior to both Lunar Flyby and 
Entry Interface (El). The TCMs are detailed in Table 2. 

The TCM performance for the three cases are presented in Table 3. 
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Cialunar EM 1 Trajectory 



Figure 2. EM-1 Trajectory 


Table 2. TCM Locations 


TCM 1 

TLI + 6 hours 

TCM 2 

TLI + 1 day 

TCM 3 

Lunar Flyby - 22 hours 

TCM 4 

Lunar Flyby - 6 hours 

TCM 5 

Lunar Flyby +18 hours 

TCM 6 

El - 2 1 hours 

TCM 7 

El - 6 hours 


11 


Table 3. 3 a TCM A V" Performance (m/s) 



One Hour Optical 

Continuous Optical 

Ground Tracking 

TCM 1 AV (m/s) 

5.14 

5.14 

5.16 

TCM 2 AV (m/s) 

0.96 

1.02 

0.19 

TCM 3 AV (m/s) 

1.55 

0.81 

0.18 

TCM 4 AV (m/s) 

0.32 

0.24 

0.31 

TCM 5 AV (m/s) 

1.53 

1.45 

1.81 

TCM 6 AV (m/s) 

3.52 

1.08 

0.66 

TCM 7 AV (m/s) 

2.80 

0.90 

0.71 

Total AU (m/s) 

15.82 

10.64 

9.03 


The Entry Flight Path Angle (FPA) delivery (at Entry Interface) for each of these cases is included 
in Table 4. 


Table 4. 3a Delivery Entry Flight Path Angle 



One Hour Optical 

Continuous Optical 

Ground Tracking 

3 a Entry Flight Path Angle (deg) 

0.252 

0.093 

0.046 


Whereas one can obtain some information regarding navigation and dispersion errors in terms 
of position and velocity, it is most helpful and illuminating to map these (instantaneous) errors 
to the entry flight path angle, because that is the quantity of ultimate relevance for being able to 
successfully complete the mission. This was done by means of state transition matrices, which were 
used to map the instantaneous navigation and dispersion errors to the time of entry interface and 
then using the partials of flight-path angle with respect to the state to complete the process. 

CONCLUDING COMMENTS 

This paper has detailed the navigation error and trajectory dispersion for the EM- 1 Free Return 
mission. It has determined that with one hour of tracking, the navigation performance is sufficient 
to meet the entry flight angle delivery error of 0.27 degrees (3 a). 
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Figure 3. Entry Flight Path Angle Navigation Error (One Hour Tracking) 
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3-or Flight Path Angle Dispersions at El (deg) 



Figure 4. Entry Flight Path Angle Trajectory Dispersion Error (One Hour Tracking) 
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3-« Navigation Errors in El Flight Path Angle (deg 



Figure 5. Entry Flight Path Angle Navigation Error (Continuous Optical Tracking) 
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Figure 6. Entry Flight Path Angle Trajectory Dispersion Error (Continuous Optical Tracking) 
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3-« Navigation Errors in El Flight Path Angle (de 



Figure 7. Entry Flight Path Angle Navigation Error (Ground Tracking Only) 
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Figure 8. Entry Flight Path Angle Trajectory Dispersion Error (Ground Tracking Only) 
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